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Shot noise in parallel wires 

Johan Lagerqvist, Yu-Chang Chen, and Massimiliano Di Ventra0| 

Department of Physics, University of California, San Diego, La Jolla, CA 92093-0319 

We report first-principles calculations of shot noise properties of parallel carbon wires in the regime 
in which the interwire distance is much smaller than the inelastic mean free path. We find that, 
with increasing interwire distance, the current approaches rapidly a value close to twice the current 
of each wire, while the Fano factor, for the same distances, is still larger than the Fano factor of a 
single wire. This enhanced Fano factor is the signature of the correlation between electron waves 
traveling along the two different wires. In addition, we find that the Fano factor is very sensitive 
to bonding between the wires, and can vary by orders of magnitudes by changing the interwire 
spacing by less than 0.5 A. While these findings confirm that shot noise is a very sensitive tool to 
probe electron transport properties in nanoscale conductors, they also suggest that a considerable 
experimental control of these structures is required to employ them in electronics. 
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INTRODUCTION 

In recent years we have witnessed an increased inter- 
est in the transport properties of molecules sandwiched 
between two bulk electrodes. H S H A, 5, d S H 
H ES m El El El E3 This interest is partly moti- 
vated by the possible use of such structures in future 
electronic applications, and partly by the desire to un- 
derstand electronic transport in atomic-scale conductors. 
In most of the experiments reported so far the number 
of molecules sandwiched between two bulk electrodes is 
unknown. US, HllllllllllllElElEi Even 
in experiments where this number is supposedly very 
small 1, 2, 3, 4, |1 H 0, |^] it is still unclear whether 
a single molecule is responsible for the reported current- 
voltage (I-V ) characteristics or if many molecules con- 
tribute. Ir^ Irrespective of this issue, future elec- 
tronic applications of molecular wires will likely require 
assembly of many components with interwire separations 
at the nanometer scale. At these distances electrons are 
less likely to experience inelastic scattering events so that 
quantum interference between electron waves travelling 
in adjacent wires dominates. We can therefore expect 
that the transport properties of an ensemble of coupled 
wires are not necessarily a simple superposition of the 
transport properties of each isolated wire. 18] Past work 
has indeed addressed this issue by studying the elec- 
tric current in parallel atomic-scale wires as a function 
of their separation. 0, 0, ^] For instance, Lang and 
Avouris have found large variations of the low-bias con- 
ductance of two carbon wires connected in parallel to 
bulk electrodes as a function of their separation. Ell 
These variations have been attributed to both a direct 
bonding between the wires and an indirect interaction 
due to the presence of the metal electrodes. El 

In this paper we want to build on the knowledge 
gained from these earlier studies and investigate a re- 
lated transport property: steady-state current fluctua- 
tions (shot noise) in parallel wires. Our motivation is 
twofold: for one, current fluctuations need to be small 



for nanoscale components to have practical applications. 
Understanding the magnitude of these fluctuations in sin- 
gle wires and to what extent interference effects between 
coupled wires modify shot noise is therefore paramount 
to progress in nanoscale electronics. Secondly, shot noise 
contains more information on the tr ansp ort properties 
of particles than the average current. |2 If In the case of 
electrons, for instance, shot noise provides information 
on their energy distribution jll, kinetics, [ll and in- 
teractions due to the Coulomb repulsion and the Pauli 
exclusion principle j!!]. Shot noise has also been found 
to depend on the atomic details of nanoscale wires like 
their length and coupling to bulk electrodes, m All these 
studies suggest that shot noise is a very sensitive tool to 
study electron transport properties in nanoscale conduc- 
tors and we hope that our work will motivate further 
experimental work in this direction. 

To the best of our knowledge very few reports have ap- 
peared in literature on the shot noise properties of paral- 
lel conductors connected to bulk electrodes. |2a, |23, |23| 
All of them employ simple models of the electronic struc- 
ture of the system and/or of the current correlations on 
each electrode. jH For instance, lannaccone et 

al. have shown that if an electron transmitted from one 
wire in a given lead does not enter the adjacent wire on 
the same lead, then the shot noise of each wire adds clas- 
sically to give a Fano factor equal to the Fano factor of 
one wire. |2y| We call this case an "uncoupled-wire" sit- 
uation. However, when the distance between the wires 
is few nanometers it is unlikely that the above case can 
occur: any phase-breaking process (whether it is due to 
inelastic scattering, quasi-elastic dephasing or electron- 
electron scattering) would not be enough to break com- 
pletely the correlation between waves traveling along the 
two distinct wires. The shot noise of the total system, 
therefore, should not be a simple classical summation of 
noise of each individual component. This is precisely the 
physical situation we want to investigate here. 

Following the work of Lang and Avouris El we study 
the transport properties of a system comprising two car- 
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bon wires of different lengths (also known as cumulene 
carbon atom chains) between two bulk electrodes (see 
schematic in Fig. The bulk electrodes are modelled 
with a semi-infinite uniform-background model (jellium 
model), m US 113 The interior electron density of the 
electrodes is taken close to the value for metallic gold 
(rg Ki 3). Note that in the work reported in Ref. ll^ the 
interior electron density of the electrodes has been cho- 
sen to have « 2. Therefore, the actual shape of the 
conductance as a function of interwire distance is differ- 
ent than the one reported in the present work. However, 
the physical trends are similar. 

The cumulene-chain system is structurally and elec- 
tronically quite simple, such that the first-principles cal- 
culations reported here can be done with a relatively high 
level of accuracy. In addition, cumulene atom chains can 
be loosely identified as the smallest "carbon nanotubes" 
possible. The results presented in this work can thus 
bear some relevance in nanotube electronics. We find, 
in agreement with previous reports, |l8l | that the current 
varies with interwire distance and saturates to a value 
close to twice the current of each wire with increasing 
wire separation. On the other hand, the Fano factor can 
vary by orders of magnitudes by changing the interwire 
spacing by less than 0.5 A. Furthermore, when the cur- 
rent is very close to the value of twice the current of each 
wire, the Fano factor is still larger than the Fano factor 
of a single wire. This enhanced Fano factor is the signa- 
ture of the correlation between electron waves traveling 
along the two different wires. At the distances consid- 
ered here, this correlation is unlikely to be destroyed by 
inelastic scattering effects, therefore this effect should be 
observed in parallel nanoscale wires. 



systems. [2^ [sO, 0| We repeat here for the reader the 
derivation of such an expression |25l | and we demonstrate 
that it reduces to Eq. when the transmission probabil- 
ities are extracted from the single-particle wavefunctions. 

We start by introducing the field operator of propa- 
gating electrons for a sample connected to a left (L) and 
right (R) electrode in terms of single-particle wavefunc- 
tions (r,K||) with energy E (atomic units are used 
throughout this paper) and component of the momentum 
parallel to the electrode surface Kii 

^^i/L^^R^ (2) 

where 

^J^(^)=.^4(«)(t)vI;^(^)(r,K||). (3) 

E 

In this expression a^^^ (t) = exp(— iwt)a^*'^'', where 
a^^^^ are the annihilation operators for electrons incident 
from the left (right) electrode, satisfying the usual anti- 
commutation relations: {a^,a'^,} = 6ijd {E — E'). We 
assume that the electrons coming from the left (right) 
electrode thermalize completely far away from the sam- 
ple, i.e., the statistics of electrons coming from left (right) 
is determined by the equilibrium Fermi-Dirac distribu- 
tion function deep into the left (right) electrode, 

i.e., 

<4ai,,>^S,,6iE-E')rE, (4) 

with i, j=R, L. 

In order to obtain the wave functions, (r, K||), we 

solve the Lippman-Schwinger equation self-consistently 



THEORETICAL APPROACH 

In the theory of mesoscopic systems, shot noise is gen- 
erally expressed in terms of the transmission probabilities 
T„ of each conducting channel n. Instead of discussing 
shot noise in a system that carries an average current /, 
it is customary to report on the ratio (called Fano fac- 
tor) between the shot noise and its Poisson (uncorrelated 
electrons) limit 2e/, where e is the electron charge. For 
a two-terminal conductor the Fano factor can then be 
expressed, in linear response, as 



F = ^T„(l-T„)/^r„ 



(1) 



vI.r(r,K,,)=<|) (r) 



I dr, J dr2G(r,ri)V^(ri,r2)*B^''^(r2,K|| 



(5) 



^0 is^K|| ^'"^ wave functions incident from the left 
(right) bare biased electrodes, and G is the corresponding 
Green's function. V is the difference between the poten- 
tial of the total system and the one of the bare biased 
electrodes. It includes the sum of the nuclear, Hartree 
and exchange-correlation potentials. 29, 30, 31] 

^From the field operator Eq. (PJ we can define the cur- 
rent operator 



However, we will not use the above expression (or 
its generalization to the non- linear case [21|). Instead, 
we will derive an expression of shot noise in terms 
of single-particle wavefunctions. |2^ This choice is a 
convenient one: we calculate the latter quantities di- 
rectly from the integral form of the Schrodinger equation 
within the density-functional theory of many-electron 
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(6) 



The ensemble average current at zero temperature is thus 



< I >= -i 



Efr 



dE / dR 



rfK|| (if ^) , 



(7) 
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where 



(8) 



with i, j=R, L. We have also assumed the right chemical 
potential to be higher than the left one. 

The shot noise spectral density is defined as the Fourier 
transform of the current autocorrelation function 



(9) 



2ttS {uo) = / dte''^' ( AI (t) AI (0) ) , 



where AI(t) is the excess current operator with respect 
to the average current defined above. Using the Bloch-De 
Dominicis theorem [32i| we can express S{uj) as 



5H= Yl fdEfh+^[l-fE) 



dRi / dK^q^^ ,, / dR2 / dK2/£^+^. (10) 



Taking the zero frequency and zero temperature limit, 
the shot noise, S, can finally be written as [l^ 



pEpR 




dE 




'J Efl 



(11) 



where is defined in Eq. (jSJ. 

In order to demonstrate that Eq. Hll() reduces to 
Eq. we consider only a single one-dimensional scatter- 
ing channel (generalization to many channels is straight- 
forward). In this case the right- and left-moving waves 
satisfy the boundary conditions (with obvious notation) 



and 
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X > oo 



(12) 



X < — oo 



(13) 
respectively. 

Introducing the above expressions into Eq. ^ and re- 
writing Eq. Hll|l as 



S{xi,X2) 



dEljij^{x,)lji'^{x2 



(14) 



we can define S^^ — S{xi — > ±00,0:2 — > ±00), with i, 
j=R, L. We thus obtain that (with Tr = Tl=T) 



:<LL 



= — r(i-r), 

TT 



(15) 



where we have implicitly assumed that Epn—EpL = eVb, 
where Vq is the external bias, and T is independent of 
energy. For T ^ 1 we obtain the Poisson limit of shot 
noise = ^T. Taking the ratio of Eq. ^ with 



leads to Eq. . Similarly we can prove that the following 
relations are valid 2l| 



S^^ = S 



RR 



LR 



(16) 



as a direct consequence of current conservation. 

As a test of Eq. pi|) we have calculated the noise for 
a single gold atom between two bulk electrodes modeled 
by jellium with — 3. The gold atom is placed at 2.2 
a.u. from each jellium edge. For a bias of O.IV we obtain 
a conductance of l.lGo (Go = 2e'^/h) and a Fano fac- 
tor of 0.14. Experimental results on gold point contacts 
indicate that to such value of conductance corresponds 
a Fano factor of 0.16 ± 0.03, |U 0] indicating a good 
agreement between theory and experiment in the present 
case. 



RESULTS AND DISCUSSION 

We are now ready to discuss the results for the par- 
allel carbon wires. We will consider two different wire 
lengths: 4-atom chains and 5-atom chains. In each wire, 
the C-C distance has been fixed at 2.5 a.u. and, due to 
the construction of the jellium model, [s^ the end car- 
bon atoms of the chains are 1.4 a.u. inside the positive 
background edge of jellium. jl^ The value of the exter- 
nal bias is y = 0.01^. The conductance in units of Go 
(top) and Fano factor (bottom) for these two systems as 
a function of the interwire distance are plotted in Figs. [21 
and 121 In order to keep the same numerical error in the 
calculation of current and noise for all interwire distances, 
and due to computational limits, we could not evaluate 
these quantities for distances larger than about 6.5 a.u. 
In these figures we report both the actual data points 
(filled squares) calculated from first principles and an in- 
terpolation in between the data points (sohd fine). We 
also indicate in the figures the "uncoupled-wire" case for 
both the current (twice the current of each independent 
wire, i.e. without any wire-wire interaction) and Fano 
factor (equal to the Fano factor of each wire - the current 
doubles as well the noise, leaving the ratio unchanged). 
It is immediately evident from both figures that both the 
current and the Fano factor change in a non-linear way 
with interwire distance. However, the variation of noise 
is even more dramatic. We discuss the current first. 

It is clear from Figs. [21 and [3 that the conductance 
variation depends on the length of the carbon wires. For 
instance, the conductance reaches a value of about twice 
the conductance of each independent wire at about 6.5 
a.u. for the 5-atom chain, while this is not the case for the 
4-atom chain. In addition, the conductance of the system 
at infinite distance (i.e. in the absence of wire- wire inter- 
actions) is larger for the 4-atom wire than for the 5-atom 
wire. This difference is related to the fact that even- 
atom-number wires have fully occupied tt states, while 
the odd-atom-number wires have a half-filled tt state at 
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the Fermi level, [l^ Esf Note, however, that the band 
width of the metal electrodes (determined by the value 
of rs) is much smaller in the present case than the one 
used in Ref.^^ This implies that the extra a and tt states 
introduced by the couphng of the wires to the electrodes 
are truly bound states in our case, thus effectively revert- 
ing the trend of the change of conductance as a function 
of wire length: in the present case, the conductance (for 
large interwire separations) is larger for even number of 
carbon atoms than for odd number of carbon atoms. Fi- 
nally, the change in conductance is related to the detailed 
bonding properties of the parallel wires. 

It is evident in Fig. |21 and Fig. O that there is a corre- 
lation between the change in conductance as a function 
of interwire distance and the corresponding change in 
Fano factor. In general, whenever the conductance in- 
creases the Fano factor decreases, as can be understood 
intuitively from Eq. However, the Fano factor is 

more affected by the interwire interaction than the con- 
ductance. We therefore discuss the bonding properties 
in connection to the Fano factor. In order to do so we 
show contour plots of the charge density. The charge den- 
sity we consider is the one due to the populated global 
current-carrying states originating from the right elec- 
trode for the system with the carbon atoms, minus the 
same quantity for the electrode-electrode system with- 
out the atoms. For the 4-carbon wire this quantity is 
plotted in Fig. 0] for a wire distance of about 2.9 a.u. 
(top) and 3.3 a.u. (bottom), i.e. for the case in which 
the noise is minimum (even lower than the value for the 
"uncouplcd-wire" case) and maximum, respectively (see 
also Fig. EJ. It is evident from Fig. |3| (top) that at 2.9 
a.u. the wires form two symmetric a bonds from each 
carbon p orbital oriented parallel to the electrode surface 
and in the plane where the centers of all C atoms lie. 
These two a bonds link the four carbon atoms between 
the jellium edges as shown in Fig. 01 (top). In this config- 
uration, electron waves travel from one electrode to the 
other as if the structure had two perfectly (degenerate) 
open channels (the value of the conductance in this case 
is also close to two, see Fig (SI . We stress however that 
this is only an analogy and does not necessarily imply 
that the number of open global channels in the structure 
is exactly two. On the other hand, at a distance of 3.3 
a.u. (bottom of Fig. there is more redistribution of 
charge from the interchain a bonds to the tt bonds along 
each wire, thus increasing backscattering and the correla- 
tion between electron waves on each side of the structure. 
This feature is also present (even if at different distances) 
for the 5-atom chain. In this case, the formation of three 
quasi-degenerate a bonds extends from about 3.5 a.u. 
to about 3.7 a.u. giving rise to a conductance of about 
3Go and a (almost) zero Fano factor in the same range 
(see Fig. . The Fano factor is thus very sensitive to 
wire-wire interactions, and in particular to the charge re- 
distribution between bonds. Since even small charge re- 



distributions can increase or reduce backscattering, the 
Fano factor can change by simply changing the interwire 
distance by less than O.SA. 

When the two wires are separated by a relatively large 
distance, the wire-wire interaction is negligible and the 
current reaches a value close to twice the value of the cur- 
rent of each wire. However, correlations between waves 
travelling in the two different wires do not necessarily 
go to zero, and this would affect the noise. This can be 
understood by looking at the contour plot of the charge 
density (as defined above) for a 5-C wire at interwire 
distance of 6.5 a.u. (Fig. EJ. At this distance the con- 
ductance is almost exactly twice the value of the conduc- 
tance of each wire (see Fig.EJ. However, the Fano factor 
is still larger than the Fano factor of each wire. This 
is because, even though there is essentially no bonding 
between the carbon atoms of the two wires (see Fig. 01, 
there is still some probability for a wave that gets out 
of a wire on one side of the junction to travel along the 
surface of the electrode and get reflected back into the 
other wire from the same side of the junction. In other 
words, the global current-carrying wavefunctions have fi- 
nite extension in the surface region between the wires. 
Indeed, for large wire separations, the lateral component 
of the scattering wavefunctions resembles closely a Bloch 
wave within few atomic units from the surface region 
in between the wires. This wavefunction extension thus 
generates extra noise compared to the "uncoupled-wire" 
case. An alternative and intuitive way of thinking about 
this point is the following: at large wire separations, the 
different paths along which an electron is carried from 
one lead to the other across either one of the two single 
wires, have energies that are quasi-degenerate. The sys- 
tem thus fluctuates fast between these quasi-degenerate 
states and additional noise is generated (in a manner sim- 
ilar to the generation of telegraph noise). Clearly, such 
correlations would be destroyed if inelastic effects were 
into play. However, as explained above, for the distances 
considered here these effects are less likely to contribute. 

We conclude by stressing that, as the present work also 
shows, shot noise is a very sensitive tool to probe electron 
transport properties. Therefore, in order to increase our 
understanding of transport in atomic-scale structures, it 
would be very desirable to have experimental results of 
noise properties of molecular wires. We hope our work 
will inspire future studies in this direction. However, 
we note that the extreme sensitivity of noise to bond- 
ing properties between parallel wires can actually consti- 
tute a limitation in electronic applications of nanoscale 
structures. 
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FIG. 1: Schematic of the system investigated. It consists of 
two parallel wires, separated by a distance L, between two 
biased electrodes. 
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FIG. 2: Conductance (upper panel), in units of quantum of 
conductance Go, and Fano factor (lower panel) as a func- 
tion of wire separation for two parallel 4-atom wires. Filled 
squares correspond to calculated values. Horizontal broken 
lines correspond to twice the conductance of a single wire (up- 
per panel) and the Fano factor of a single wire (lower panel) . 
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FIG. 3: Conductance (upper panel), in units of quantum of 
conductance Go, and Fano factor (lower panel) as a func- 
tion of wire separation for two parallel 5-atoni wires. Filled 
squares correspond to calculated values. Horizontal broken 
lines correspond to twice the conductance of a single wire (up- 
per panel) and the Fano factor of a single wire (lower panel). 
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FIG. 4: Contour plot of the charge density for the 4-atom wire 
system at two different interwire distances. See text for the 
exphcit definition of this density. Vertical black lines corre- 
spond to the edges of the jellium model and circles correspond 
to atomic positions. 
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FIG. 5: Contour plot of the charge density for the 5-atom 
wire system at interwire distance of 6.5 a.u. See text for the 
explicit definition of this density. Vertical black lines corre- 
spond to the edges of the jellium model and circles correspond 
to atomic positions. 



